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Abstract 

In this letter we present our results for the four-loop beta-function of the strong coupling in the Standard 
Model of fundamental interactions. The expression is obtained from gluon self-energy diagrams in the 
background field gauge without application of special infra-red rearrangement tricks. We take top-Yukawa 
and self-Higgs interactions into account, but neglect electroweak gauge couplings. Ambiguities due to 75 
treatment are discussed and a particular “reading” prescription for odd Dirac traces is advocated. 
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Most of viable models describing Nature at high energies are based on gauge symmetries. Quantum 
chromodynamics (QCD) is a gauge theory of strong interactions and it is important to study its strength 
described by the coupling a s both in the low- and high-energy limit. At low energies the QCD interactions 
play a dominant role in binding quarks and gluons together into nucleons. At larger scales the coupling a s 
decreases m due to non-Abelian nature of the underlying gauge theory. Nevertheless, the precise value of 
a s is of paramount importance both for modern and future colliders and for theoretical studies of physics 
going beyond the Standard Model (SM) of fundamental interactions. 

The renormalization group equations (RGE) relate couplings at different scales. By solving them one 
can not only confront measurements carried out at different energies with theory but also study asymptotic 
behavior of the latter at scales inaccessible in current and even future experiments. 

The progress in calculation of beta-functions — key RGE quantities — is tightly connected to the 
introduction of dimensional regularization [3j and (modified) minimal subtraction (or MS) scheme. The 
former does not break gauge symmetry in d = 4 — 2e dimensions and the advantage of the latter lies in 
the fact that beta-functions are extracted only from ultraviolet (UV) asymptotics of Feynman integrals. 
This fact allows one to drastically simplify a calculation by modifying the infra-red (IR) structure of the 
considered integrals by means of the so-called infra-red rearrangement (IRR) procedure [4] . 

Pure QCD results for the strong coupling beta-function are known for quite a long time up to four 
loops mmmm. The four-loop results are obtained with the help of IRR procedure leading to four-loop 
vacuum diagrams with all lines having the same mass. The beta-function is extracted from the quark-gluon 
uni and ghost-gluon m vertex renormalization constants and the corresponding wave function renor¬ 
malization constants. Independently, the ghost field renormalization constant and the ghost-gluon vertex 
renormalization constant were calculated in Ref. |12j . The quark field renormalization constant j!3] was 
found by a different method bringing the problem to the calculation of three-loop massless propagator-type 
diagrams. 

It is obvious that in a precise study of QCD processes at high energies one should consider other SM 
interactions and their effect on the running of a s . Recently, the full set of three-loop beta-functions for all 
SM parameters, including the strong coupling, was found in a series of papers |14H20j . 
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In this letter we present our result for the dominant four-loop contribution to the beta-function of the 
strong coupling in the SM. In our calculation we neglect the electroweak gauge interactions, but keep top- 
Yukawa and Higgs self-interactions along with the well-known QCD corrections. The calculation is carried 
out in the background-field gauge (BFG) [22 [22] ■ The advantage of BFG lies in the QED-like relation 
between the gauge coupling renormalization constant Z as and that of the background gluon field ^3 : 

00 X7(i) 

z aB = i/Zy 3 , z = l + 'Yl,—r- • (!) 

i=1 

This allows us to obtain the final result by considering massless propagator-type integrals. For the present 
calculation we evaluate the required four-loop integrals contributing to the two-point Green function for the 
background gluon field and exploit multiplicative renormalizability of theory as in Refs. pencil]. We only 
need to calculate bare Green-functions up to the Z-loop level and re-express all the model parameters in 
terms of the renormalized ones in the MS-scheme: 
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where a collectively denotes the SM couplings in the considered limit, i.e., strong - g s , top-Yukawa - y t , and 
that of Higgs self-interaction - A, together with SU(3) gauge-fixing parameter £: 

(167r 2 )a = {^,y t 2 ,A, (16 tt 2 )C} • (3) 

As it is clear from ([ 2 ]), due to finiteness of r ren , RHS should be finite too and at each order of perturbation 
theory we have an equation on Zr terms. Given the knowledge of all three-loop renormalization constants, 
it is possible to extract the four-loop contribution to the renormalization constant Zy 3 for the background 
gluon field V3. 

For diagram generation the package DIANA [25], which internally uses QGRAF [25], was utilized. After 
some color |27] and Dirac algebra all the generated two-point functions were mapped onto three auxiliary 
topologies, each containing 11 propagators and 3 irreducible numerators. Before the actual integration-by¬ 
parts (IBP) [28] reduction, a set of reduction rules was prepared by means of LiteRed{29] package. After 
that the prepared set was passed to the C++ version of FIRE package [30] , which allows one to reduce the 
obtained four-loop integrals in parallel on a multithread machine. The IBP reduction leads to a small set of 
master integrals. The expressions for the latter are known in analytical form up to the finite parts [5T| . The 
master integrals were also cross-checked numerically in Ref. 1321 and some additional terms in e-expansion 
were found in Ref. [33]. We also perform some independent evaluation of several three-loop integrals. The 
simplest four-loop integrals with three-loop insertions in a one-loop integral were checked by means of the 
FORM-based package MINCER [MH55] . 

It is also worth mentioning that as an independent cross-check of our setup, we prepared a simple 
QCD model with additional fermion in the adjoint representation of SU(3) color group (“gluino”). The 
beta-function for such a model at the three-loop order can be predicted by means of proper color factor 
substitutions [36] , At four loops similar procedure ceases to be sufficient and some additional information 
from the direct calculations m of QED-type diagrams is required to predict the beta-function. We compared 
the predicted results [35] with that calculated by means of the above-mentioned setup and found a perfect 
agreement. 

After these kind of tests we address the issue of finding the SM four-loop contribution to the gauge 
coupling beta-function. The calculation of bare Green-functions requires evaluation of traces over Dirac 
matrices in d ^ 4 dimensions and an additional complication, comparing to QCD, arises when traces 
involving y 5 matrices are present (see, e.g., Ref. 59] ). 

One needs to be careful in maintaining 75 anticommutativity and strict four-dimensional relation 


tr (7 / V7 P 7 a 7 5 ) = -He^, 


(4) 
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involving totally antisymmetric tensor with e 0123 = 1. At three loops m ns it was proven by direct 
calculations that it is possible to use senri-naive approach and utilize both { 7 ^, 75 } = 0 and Q without 
paying much attention to apparent non-cyclicity of the trace operation in d / 4. The formal e-tensor 
originating from the trace 0 can only give a non-trivial contribution if two such traces are present and the 
anti-symmetric tensors are contracted by means of 

e^°e afhS = W- (5) 

Strictly speaking, the Kronecker delta-symbols in ([ 5 ]) should be considered as four-dimensional objects and 
the contraction with the remaining part of a diagram should be carried out after subtraction of infinities via 
A-operation |41j . However, it is not convenient in a massive calculation involving thousands of Feynman 
diagrams and it is tempting to use ^-dimensional satisfying <5^ = c? in the bare d = 4 — 2e theory. It 
is easy to convince oneself that the cyclic property of traces gives rise to an ambiguity 0(e), which can 
play a role in determining RG coefficients. However, it turns out that at the three-loop level a non-trivial 
contribution originating from the contraction ([ 5 ]) appears for the first time in the Yukawa coupling beta- 
function [40| . Both gauge-coupling and higgs self-coupling turn out to be free from this kind of contributions 
due to gauge-anomaly cancellation conditions fulfilled with the SM. 

In our calculation we tried to employ the above-mentioned semi-naive approach to study the diagrams 
giving rise to a non-zero terms due to ([ 5 ]). A typical diagram is shown in Fig. [I] By counting coupling 
constants and performing color algebra it is easy to convince oneself that 75 affects only a 2 a|T|. contribution 
with Tp = 1/2. We have 24 planar diagrams of this type and 48 non-planar graphs, which can be obtained by 
permutation of internal lines connecting two fermion traces. It contains box-type closed sub-loops involving 
two gauge bosons and two scalars as external legs [^] Direct evaluation of the diagram shows that, indeed, 
the resulting expression is free from higher poles in e. 

An additional argument to the fact that the four-loop higher poles are not affected by different 75 
prescriptions comes from an observation that they can be found in advance from the known three-loop 
results fl4f[20l and the so-called pole equations [3j. Dangerous contributions due to ([5]) to the four-loop 
higher poles can only appear if three-loop Yukawa coupling beta-function p a± is involved. However, it is 
easy to prove that it is not the case for f} ag • Since one-loop gauge-coupling beta-functions does not depend 
on other couplings, only one- and two-loop terms in /3 at contribute to the four-loop pole equations for a s 
and we expect no dangerous high-order poles in Z ag . It is worth stressing that the argument can also be 
applied to all gauge-coupling beta-functions in the full SM provided that all gauge anomalies are canceled. 

This kind of reasoning lead us to a premature conclusion that the semi-naive approach is sufficient to get 
the correct answer for the strong coupling beta-function at four loops. The question whether the absence of 
higher poles in the considered diagrams is sufficient for them to be unambiguous was initially left without 
consideration. 

However, a more careful study, triggered by the appearance of similar result in Ref. [42], of all the 
above-mentioned 72 diagrams have shown that our initial treatment leads to results inconsistent with formal 
charge-conjugation symmetry. For example, the left fermion loop in FigjT] gives 

tr (Pi7pP2P375P47/x) = tr (Pl7pP2P375P47p) T , ( 6 ) 

where p, p are Lorentz indices of gluons connected to the loop, Pi correspond to momenta of fermion prop¬ 
agators, and 75 comes from the vertex with neutral would-be goldstone boson (%). Assuming standard 
relations C^C ” 1 = — 7 J and C75C -1 = 7 J due to charge conjugation transformation C , one can show 
that the right-hand side (RHS) of Eq. § can be rewritten in the form 

tr (7 p(-P4)75(-P3)(-P2)7p(-Pi)) > (7) 

which corresponds to a diagram with flipped fermion flow of the considered closed chain. As a consequence, 
these two diagrams should be equal. However, an inconsistent treatment of 75 could spoil this C-symmetry. 


2 Similar sub-diagrams of three-loop Yukawa vertex were discussed in papers [16:. 40l in the same context. 
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As it is knowJJS], the cyclic property of the Dirac trace should be abandoned if one tries to utilize both 75 
anticommutativity and trace relation Q. Due to this, the “reading” point [43], which defines the position 
of 75 in an “odd” trace, should be chosen consistently for diagrams related by C-symmetry. In our initial 
treatment the choice was made by the utilized diagram generation code DIANA, which starts writing closed 
fermion loops from a propagator with certain momentum, irrespectively of its fermion flow, e.g., from (— pi) 
for Eq. ([7]). As it will be shown latter, this prescription leads to an inconsistent result. 

In order to shed light on possible solutions to this issue, i.e., we study different “reading” points [13] for 
the involved Dirac traces. We “cut” a trace at a certain point, i.e., start writing a fermion string from either 
propagator or vertex, and { 75 , 7 M } = 0 , 7 ! = 1 properties to anticommute 75 to the rightmost position of 
the corresponding chain. For example, Eq. (| 6 | corresponds to cutting the trace at the external gauge vertex, 
while the flipped version ([TJ) is cut at the propagator entering the vertex. The resulting traces involve more 
than four 7 -matrices and direct application of the trace condition Q is not possible. We made use of the 
algorithm given in the FORM [44] manual to reduce these traces to Eq. Q. 

Simple “scan” over different reading points indicates that, indeed, there is an ambiguity in the resulting 
expression depending on the position, at which the traces are “cut”. We can distinguish three situations: 
A) both traces start (or encQ at external vertices; B) only one trace starts(or ends) at an external vertex; 
C) both traces start(or end) at some internal vertex. The ambiguous divergent contribution to the gluon 
self-energy from a diagram involving odd traces with 75 can be parameterized in the following way 

n 2 n 2 T 2 

a ^-^(X 1 +X 2 C 3 ). ( 8 ) 

e 

For the case A, the planar graphs give X\ = 1/6, X 2 = 0, while non-planar ones lead to X\ = —1/18 and 
X 2 = 1/6. Summing contributions from 24 planar and 48 non-planar diagrams we obtain 


a 2 a 2 T 2 


€ 



(9) 


For the case B we have found that the corresponding coefficients are multiplied by a factor of 2, while 
the prescription C give rise to a factor of 3. 

It is worth mentioning that we also tried to utilize Larin-like [45] prescription for dealing with 75 , i.e., 
the substitution 

7m 75 -t ^nupcrYYl' 7 , (10) 

which, in a cyclic trace with anticommuting 75 , can be interpreted as a reading prescription, which corre¬ 
sponds to an average of two adjacent cut-points 


tr (■ • • 7 ^ 75 ) - [tr (... 7 ^ 75 ) - tr ( 7m ... 75 )] . 


( 11 ) 


This latter fact was confirmed by direct calculations. There is a subtlety with (10) coming from the fact that 


7 m may not only come from a gauge vertex, but from a propagator, connecting external and internal vertices. 
For example, if one terminates the left fermion chain in Fig. [T] by the external vertex, but the right one — 
by the Yukawa vertex involving, e.g., a neutral Goldsone boson \ (marked by a blue circle), the outcome 
will be the mean of the results corresponding to the cases A and B, i.e., X\ = 1/2 • 1/6 • {\a + 2 b) = 1/4. 

Similar situation appears if we do not move 75 in (a part of) the trace with single 75 , but just use fl0| 
without anticommuting it to some cut-point. This approach turns out to be equivalent to the one, when 
the cut-point is actually fixed to be the vertex, in which 75 appears. Again, referring to Fig. [l]with internal 
X, application of (101 from the very beginning gives rise to X\ = 1/2 • 1/6 • (2 b + 3 c) = 5/12. The only 
complication in the charged (jX goldstones case comes from the necessity of additional averaging over four 


3 We can also consider situations when traces are terminated at fermion propagators, but it is easy to convince oneself that 

these cases are also included in our consideration. 
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different ways to choose the position of single 75 in both traces, since 1 ± 75 enters each vertex with cjr 1 . 
For example, 

Al = 22 ' 22 ' 6 ’ ^’ 4 + 3c ’ 3 ’ 4 ^ = 24 ’ ( 12 ) 

However, it is easy to understand that such a prescription will give rise to a result incompatible with the 
above-mentioned C-symmetry, since flipping the right loop fermion flow results in 

Xl = 6 ' (3c ' 4 ' 4) = 2' (13) 


Nevertheless, it is worth pointing that the issue arises only due to identical handling of both diagrams. In 
general, one is able to use different reading prescriptions for different diagrams and by a proper choice it is 
possible to maintain the C-symmetry. We can parametrize our ignorance of the correct reading prescription 
by means of a function 


R(x, y, z) = 


1 A ' x + 2 B ' y + 3c ■ z 
x + y + z 


1 < R(x, y, z) < 3, 


which corresponds to an average over x + y + z reading points, x of which are 
z - of type C. 


x,y,z£ N, (14) 

of type A, y - of type B, and 



Figure 1: A typical diagram with two fermion traces giving rise to a non-zero contribution to fta 8 at four loops due to the 
appearance of 75 in Yukawa-type vertices of the would-be Goldstone bosons y, cA. Fermion line momenta p, are indicated. 
The momentum flow in the loop coincides with the fermion flow. The latter is denoted by arrows. 


In spite of the observed ambiguity, in what follows we provide some arguments leading us to a conclusion 
that the prescription C with R = 3 is preferred among the others. 

First of all, as it is noted in Ref. [43], reading points at external vertices could spoil gauge invariance 
of the final result of our two-point function. Indeed, we have checked that finite parts of the (sum of) 
corresponding diagrams are not transverse if both Dirac traces are cut at the external vertices. It turns 
out that the prescription B also yields zero upon multiplication by the product of external momenta q^q^- 
Nevertheless, if we want to avoid non-symmetric treatment of the external vertices, we are left with the only 
reasonable prescription, which is C in our case. 

One more argument comes from utilization of Fierz identities, which convert two closed fermion chains 
into a single loop. A well-known subtlety lies in the fact that the relations in non-integer dimensions involve 
infinitely many terms (see, e.g., Ref. j46j 1. In our calculation we have tried to incorporate the following 
identity 


/ y Taia 4 Ta3a2 
p.=l 


= 2- d/2 E(-l) 


1=0 


10 + 1 ) (d — 21) 

Z! 


x aia 2 x a 3 




(15) 


in which totally anti-symmetric combinations 






= )£(-»' 


ip.i 1 ip.i 2 ■■■'y p.i 


(16) 


are introduced. In our case the two contracted indices in RHS of (151 correspond to the gluon propagator 


connecting two fermion traces in the Feynman gauge. We have checked that up to terms l = 4 the obtained 
result coincide with the one corresponding to R = 3. 
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Finally, we re-calculated the problematic diagrams by applying an infra-red rearrangement technique 
0103. i- e -> by transforming the integrals into fully-massive bubbles. An important difference from Ref. 142; 
lies in the fact that we perform Dirac traces applying self-consistent BMHV-algebra [a Eg after (tensor) 
integrals are evaluated in d-dimension. We have used spinney package j49| to keep track of the dimension 
of the involved Lorentz indices and cross-checked that the final expression does not depend on the cut point. 

After providing this kind of arguments, let us now proceed with our result. We define the four-loop 
beta-function for a s as 


d a s 

d log /x 2 


3 

Pa a a s = —a s Pih l+2 

i =o 


(17) 


where we use h to count the powers of coupling constants given in <§. Our final result for p 3 can be written 
in terms of SU(3) casimirs and the number of SM generations - n G : 


@3 = P3 CD (. n f = 2n G ) + a^at 
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00 
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+ TpCp (117 - 144C 3 ) + 222 T F C A 


V 3 
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(18) 


where pure n/-flavour QCD contribution /1^ CD is available from Refs. llQlITT], in which the same notation for 
beta-function and coupling constants was used. The result is free from gauge-fixing parameter dependence. 
This serves as a welcome check since at the intermediate steps one has to take into account that the gauge- 
fixing parameter £ is renormalized in the same way as quantum gluon field. In (18) we emphasize terms 
O {a 2 s a 2 J V 2 ), which are affected by contributions from diagrams similar to that given in Fig. [ljand multiply 
the corresponding terms by the factor R from Eq. (14). 



Figure 2: Relative size of the calculated four-loop contributions and the pure QCD five-loop term 0(7/.)) (found recently in 
Ref. [50] and indicated by arrow) with respect to their sum. Nontrivial part due to 75 is also indicated. The numbers in the 
inset should be multiplied by 10. Negligible correction due to a s a,t<i\ is not shown. 

Now we turn to the numerical impact of the addition of new terms (see Fig. [ 2 ]). Using as input the 
values of strong, Yukawa and Higgs self-coupling from [51] at the scale fj, = M t and also evolving the same 
input upto the Planck scale, we see that at both scales the leading pure QCD contribution O(a^) accounts 
~ 93% of the total four-loop result while the dominant 0(a? s a t ) terms due to top-Yukawa only account 
~ 10%. We observe that the 0(a 2 s al) contribution has an opposite sign and is about ~ 5% of f) 3 but terms 
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0(a s af) have the same sign and contribute to the order of ~ 1%. The other contributions are less than 
1% at both the considered scales. Terms due to non-trivial 75 treatment lead to about —0.5% reduction in 
fa. Nevertheless, the corresponding contribution is much larger than pure QCD five-loop corrections, found 
recently in Ref. [5T| and partially confirmed in Ref. [52]. 

In conclusion, we calculated four-loop beta-function for the SM strong coupling in the limit of vanishing 
electroweak gauge-couplings. By means of the background-field gauge formalism the result was obtained 
solely by considering four-loop massless propagator-type diagrams. We studied the 75 -issue and found that 
the senri-naive treatment is not sufficient to obtain unambiguous expression. In our approach the ambiguity 
manifests itself in the dependence of the result on the “reading” prescription, chosen to evaluate traces 
with odd number of 75 , and was parametrized by the factor 1 < R < 3 in our final expression for fa (18). 
Nevertheless, we provide several arguments supporting the prescription with R = 3, for which both traces 
are “cut” at some internal vertex. 

The obtained result is a necessary step towards a four-loop RG analysis of the SM. However, one should 
keep in mind that proper inclusion of three-loop non-QCD decoupling corrections is required to evaluate 
a s consistently at the electroweak scale in the full SM (for pure QCD results up to the four-loop level see 
Refs. [55H55] and for the leading two-loop electroweak corrections see Ref. [55]). In addition, the results can 
be applied to the study of the so-called Weyl consistency conditions foTboiTl within the SM, thus, extending 
the analysis of Ref. [60] to the four-loop level. 
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